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SUMMARY 


A  series  of  experiments  have  been  conducted  to  supplement  our  current 
data  bank  of  helmet  loading  configurations  (then  at  fifteen  different  helmet 
weight  and  center-of-gravity  combinations).  The  neck  muscles  were 
dynamically  and  statically  loaded  by  systematic  variation  of  nine  additional 
headgear  configurations  consisting  of  three  different  combinations  of 
centers-of-gravity  (right-forward-low,  left-lateral-low  and  right-af tward- 
low)  and  three  different  weights  (1.45  kg.,  2.27  kg.  and  4.09  kg.).  Six 
subjects  would  rotate  their  heads  later*' 1  ly  (from  side-to-side)  for  30 
minutes  with  each  of  the  headgear  loading  combinations.  Immediately 
thereafter,  the  subject  would  position  his  head  in  an  isometric  head 
dynamometer  and  exert  a  sustained  right  lateral  neck  contraction  or  forward 
neck  contraction  at  70%  of  his  maximum  strength,  during  which  endurance  time 
(to  fatigue)  was  recorded,  the  EMG  over  the  right  sternocleidomastoid 
muscle,  over  the  posterior  trapezius/splenius  muscles,  and  the  systolic  and 
diastolic  blood  pressure  and  heart  rate  were  continuously  recorded.  Of  the 
twenty-four  resultant  combinations  ot  helmet  weights  and  centers-of-gravity: 
(1)  18  combinations  were  used  as  the  boundary  conditions  for  an  empirical 
mathematical  model  to  predict  both  forward  and  lateral  neck  muscle  endurance 
for  any  weight-C.G.  combination  within  the  boundary  conditions,  (2)  3 
combinations  were  used  to  test  the  assumption  of  insensitivity  to  vertical 
loading;  and  (3)  3  combinations  were  used  to  test  the  assumption  of 
bilateral  symmetric  response.  Finally,  a  statistical  analysis  program  was 
supplied  to  perform  paired  T-tests  for  comparison  of  different  headgear 
loading  configurations.  The  results  indicate  that  the  mathematical  model 
makes  reasonable  predictions  within  the  boundary  conditions.  Input  data 
outside  the  boundary  conditions  is  rejected.  The  assumption  of 
insensitivity  to  vertical  loading  is  demonstrated.  The  assumption  of 
bilateral  symmetric  response  was  confirmed  for  the  1.45  kg.  and  2.27  kg. 
helmet  loads.  However,  this  assumption  was  not  confirmed  for  the  4.09  kg. 
helmet  load.  It  is  concluded  that  the  computer  model  is  valid  for  midline, 
vertical  and  lateral  headgear  loading  within  the  boundary  conditions 
specified.  (This  work  was  supported  by  U.S.  Army  contract 
DAMD17-80-C-0089. ) 
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INTRODUCTION 


This  section  will  outline  the  objectives  of  the  research  and  its 
significance  to  the  U.S.  Army.  The  previous  work  by  the  investigators  will 
be  reviewed. 

OBJECTIVES 

The  following  objectives  were  addressed  during  the  course  of  this  study: 

A.  Supplement  our  current  data  bank  of  helmet  loading  configurations 
(then  at  15  different  helmet  weight  and  center-of-gravity 
combinations)  with  additional  experiments  (9  additional  helmet 
weight  and  center-of-gravity  combinations)  in  order  to  define  the 
necessary  boundary  conditions  for  a  realistic  mathematical  model; 

B.  Develop  an  appropriate  empirical  mathematical  model  (using  the 
technique  of  piece-wise  linear  analysis)  to  predict  both  forward 
and  lateral  neck  muscle  endurance  for  any  weight-C.G.  configuration 
within  the  boundary  conditions; 

C.  Deliver  to  the  U.S.  Army,  the  results  of  this  research  in 
•'software"  package  form  that  is  (1)  conveniently  formatted,  (2) 
easily  accessible,  and  (3)  readily  interpretable. 

SIGNIFICANCE 

Our  results  are  directly  applicable  to  several  objectives  of  the  U.S. 
Army  Aeromedical  Research  Laboratory  (USAARL) .  It  will  provide  for 
systematic  understanding  of  how  different  headgear  designs  load  the  cervical 
muscles,  affect  the  fatigue  end-point,  and  therefore  the  subject’s  ability 
to  tolerate  various  headgear  loading.  Consequently,  it  will  enable  better 
(or  optimal)  helmet  design  wnicli  wnl  mininuze  both  muscle  loading  and 
fatigue  in  the  various  driving  or  flying  environments,  and  therefore 
maximize  operational  endurance  time. 

Pilots  are  currently  being  asked  tc  wear  and  use  additional  headgear. 

For  example,  night  vision  goggles  are  now  being  worn  in  combination  with  a 
helmet.  An  objective  evaluation  of  the  cumulative  effects  of  various  helmet 
weights  and  center-of-gravity  combinations  on  necx  muscle  lei.oiou  ai.d 
fatigue  is  now  needed  to  establish  the  optimal  "trade  off"  between  headgear 
requirements  and  physiological  capabilities. 

Finally,  the  U.S.  Army  continues  to  design  and  evaluate  new  helmets  for 
crew  members.  Impact  protection,  noise  protection  and  visual  protection 
(among  other  parameters)  are  all  capable  of  objective  quantitative 
evaluation  with  respect  to  helmet  design.  Successful  application  of  this 
experimental  project  will  yield  a  system  of  equations  that  will  allow  the 
designer  to  input  important  helmet  design  parameters  (i.e.,  weight  and 
center-of-gravity),  and  the  equations  will  output  isometric  endurance  time 
for  the  neck  muscles  in  the  forward  and  lateral  contraction  mode. 

BACKGROUND 

In  our  laboratory,  support  from  Army  contract  DAMD17- 80-C- 0089,  began  in 
June,  1980.  The  Biomechanics  Laboratory  at  Wright  State  University  has 
subsequently  developed  a  unique  capability  to  evaluate  neck  muscle  endurance 
and  fatigue  as  a  function  of  helmet  loading.  In  pursuing  our  studies  for 
the  Army,  two  unique  pieces  of  equipment  were  ueveloped.  First,  an 
isometric  helmet  dynamometer  was  designed.  This  allows  us  to  quantitate  the 


strength  of  neck  muscle  contractions  in  the  forward,  backward  and  both 
lateral  modes  (Phillips  and  Petrofsky  1981b,  1981c).  Furthermore,  it  allows 
subjects  to  hold  a  "target"  tension  less  than  their  maximal  contraction 
tension  while  simultaneously  measuring  (1)  endurance  time,  (2)  necl:  muscle 
EMG,  and  (3)  cardiovascular  parameters.  Second,  a  variable  weight  and 
variable  C.G.  helmet  simulator  was  developed.  This  helmet  simulator  has 
been  calibrated  for  five  weights  (1.4  to  4.1  kgs.)  and  twenty  one 
center-of-gravity  locations  (in  the  X,  Y  and  Z  plane).  Furthermore,  the 
helmet  simulator  can  be  recalibrated  for  any  desired  weight  and  C.G. 
combination  within  its  design  limits.  The  helmet  simulator  was  validated 
during  the  first  year  of  our  army  contract  (Phillips  and  Petrofsky  1982b), 
and  during  our  second  year  (commencing  June,  1981),  allowed  us  to  observe 
the  effects  of  fifteen  different  weight  and  C.G.  combinations  on  neck  muscle 
endurance  and  fatigue. 

Between  June,  1980  and  May,  1981  a  number  of  important  parameters  were 
evaluated.  For  the  first  time,  our  group  determined  the  basic  strength- 
endurance  curves  for  neck  muscle  in  forward,  right-lateral,  and  backward 
contraction  modes  and  compared  them  to  the  hand-grip  muscles  in  the  same 
subjects  (Phillips  and  Petrofsky  1981a;  Petrofsky  and  Phillips  1982).  The 
research  program  then  investigated  the  effects  of  no  helmet  (CONTROL),  a 
standard  SPH-4  helmet  (HELMET),  and  the  above  helmet  combined  with  Night- 
Vision-Goggles  (NVG)  on  the  cardiovascular  responses  (Phillips  and  Petrofsky 
1982c),  the  neck  muscle  isometric  endurance  time  and  EMG  response  (Phillips 
and  Petrofsky  1982a,  1983a).  Finally,  the  variable  center-of-gravity  and 
variable  weight  helmet  simulator  (VCGW)  was  validated  against  the  SPH-4 
helmet  and  Night-Vision-Goggle  combination  (H/NVG)  with  respect  to  neck 
muscle  endurance  time  and  EMG  response  (Phillips  and  Petrofsky  198?h),  The 
final  contract  report  has  not  yet  been  released  for  general  distribution 
pending  some  technical  and  format  revisions  (Phillips  and  Petrofsky  1981d). 

Between  June  1901  and  May  1982,  the  isometric  endurance  time,  EMG 
response,  cardiovascular  response  (Phillips  and  Petrofsky  1983c,  1984b)  ar.d 
strength-recovery  response  were  evaluated  as  a  function  of  fifteen  different 
helmet  stimulator  configurations.  These  represent  a  combination  of  three 
different  helmet  Weights  (i.e.,  1.4,  2.0,  4.1  kg .  >  and  five  different 
centers-of-gravity  (i.e.,  centav-low  (CL],  center-high  (CH) ,  forward-low 
[FL],  rearward-low  (RL)  and  right-lateral-low  [RLL] ,  as  shown  graphically  in 
Figure  1.  After  analysis  of  these  results,  it  was  concluded  that  there  were 
optimal  C.G.  locations  for  1.4  and  2.3  kg.  ( furwai d- low  and  tight- lateral 
low)  and  for  4.1  kg.  (af tward- l ow)  (Phillips  and  Petrofsky  1983b).  The 
final  contract  report  for  this  period  was  submitted  in  November,  1982,  and 
is  still  undergoing  technical  and  editorial  review  by  the  U.S.  Army 
(Phillips  and  Petrofsky  1982d). 

It  can  be  appreciated  that  any  helmet  system  (whether  current  or 
projected)  will  have  its  own  unique  weight  and  C.G.  location  which  will 
rarely  fit  one  of  the  15  combinations  studied  during  the  second  year’s  work. 

Furthermore,  it  would  be  highly  impractical  to  test  the  large  number  of 
configurations  required  so  that  any  present  or  projected  helmet  system  could 
be  closely  approximated. 

Referring  again  to  Fig.  1,  a  coi.*t..t  ional  three-dimensional  reference 
system  is  shown  and  those  configurations  noted  which  were  determined  at  the 
conclusion  of  second-year  study  (O)  and  during  the  first  half  of  the  present 
study  (•) ,  see  Methods  and  Materials  section.  The  arrow  notation  refers 
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to  positive  movements  about  each  axis,  and  numbers  beside  each  C.G. 

Co: figuration  are  weights  (in  pounds). 

The  additional  data  points  (•)  provided  the  additional  boundary 
conditions  necessary  to  apply  piece-wj se  linear  analysis.  The  resultant 
equations  are  capable  of  predicting  the  useful  operational  endurance  time  as 
a  function  of  any  helmet  weight  and  center-of-gravity  combination  within  the 
boundary  conditions  (Phillips  and  Petrofsky  1984a,  1986). 

METHODS  ADD  MATERIALS 

This  section  describes  how  the  study  was  conducted,  including  use  of 
subjects,  lhis  section  also  includes  a  brief  description  of  materials  and 
apparatus  used  in  the  study.  The  fifteen  headgear  combinations  reported  for 
the  second  year  study  are  combined  with  the  nine  headgear  combinations 
evaluated  in  this  final  year  study  (a  total  of  twenty-four  headgear 
combinations).  This  was  done  since  the  objective  of  the  final  year  study 
was  a  mathematical  model  which  required  all  twenty-four  headgear 
configurations  in  order  to  derive  the  model. 

SUBJECTS 

Six  subjects  were  used  in  these  experiments.  The  subjects  were  male 
volunteer  university  students  whose  ages,  heights,  neck  sizes,  and  weights 
are  listed  in  Table  1.  All  subjects  were  informed  of  all  experimental 
procedures  and  were  medically  examined  including  a  thorough  history  and  a 
complete  physical  exam.  All  procedures  were  fully  approved  by  the  committee 
on  human  experimentation. 

TRAINING 

All  subjects  were  first  trained  to  produce  a  maximum  voluntary  effort 
and  to  sustain  that  effort  to  fatigue  at  the  tension  used  in  the  study  and 
with  the  various  muscle  groups  examined  here.  Isometric  training  consisted 
of  a  series  of  brief  (<  3  sec.)  maximal  voluntary  contractions  (MVC)  with  an 
intercontraction  interval  of  3  minutes.  These  were  followed  by  a  fatiguing 
isometric  contraction.  The  tension  exerted  during  the  fatiguing  contraction 
was  set  at  70X  of  the  MVC.  On  any  one  day,  only  one  direction  of 
contraction  was  performed  and  all  fatiguing  contractions  were  held  at  the 
same  percentage  of  isometric  tension.  This  procedure  was  repeated  on 
Monday ,  Wednesday,  and  Friday  of  successive  weeks  until,  for  any  one  muscle 
group  (direction),  the  coefficient  of  variation  (standard  deviation  divided 
by  the  mean)  of  endurance  from  day  to  day  was  reduced  to  less  than  5%.  In 
practice,  the  coefficient  of  variation  of  strength  in  these  trained  subjects 
was  less  than  3%  from  day  to  day  by  the  end  of  the  training  period. 

Training  was  conducted  at  70%  MVC  and  with  both  muscle  groups  examined 
here.  For  most  subjects,  training  for  any  one  muscle  group  averaged  about  3 
weeks . 

ISOMETRIC  HEAD  DYNAMOMETER 

A  helmet  dynamometer  has  been  developed  which  can  be  used  to  measure  the 
strength  and  endurance  of  neck  muscles  in  man  in  either  one  of  four 
directions  (forward  flexion,  backward  extension,  right  and  left  lateral 
flexion).  The  dynamometer  is  based  around  the  army  SPH-4  type  helmet,  but 
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is  easily  adaptable  to  other  types  of  military  helmets  as  well.  The 
dynamometer  maV.es  It  possible  to  evaluate  the  effect  of  various  types  of 
dynamic  activities  ar.d  other  flight  activities  on  neck  muscle  strength  and 
neck  muscle  endurance.  It  is,  therefore,  a  useful  tool  in  the  study  of 
military  helmet  design  and  evaluation  of  the  stress  induced  by  flight 
maneuvers.  The  isometric  head  dynamometer  has  been  described  in  detail  by 
Petrofsky  and  Phillips  (1982)  and  is  shown  in  Figures  2,  3  and  4. 

HELMET  SIMULATOR 

The  systematic  assessment  of  significant  helmet  design  parameters 
employed  a  helmet  simulator  In  which  both  the  weight  and  center-of-gravity 
were  methodically  and  controllably  altered.  Such  a  helmet  was  developed  by 
Simula,  Inc.,  under  subcontract  to  Wright  State  University. 

The  helmet  simulator  consists  of  two  weight  concealment  boxes  attached 
to  opposite  sides  of  a  support  ring  (headring)  which  in  turn  is  si  'ported 
upon  the  wearer's  head  by  a  suspension  system  taken  from  an  SFH-4  helmet. 

The  weight  and  C.G.  can  be  altered  by  positioning  variable  weights  within 
the  concealment  boxes.  Fabric  covers  over  the  boxes  prevent  the  test 
subjects  from  obtaining  visual  clues  as  to  the  C.G.  location. 

The  minimum  weight  of  the  helmet  simulator,  without  any  variable  weights 
in  the  boxes,  is  2.5  lb.,  slightly  less  than  the  weight  of  most  quality 
crash  helmets  made  by  reputable  manufacturers.  The  addition  of  variable 
weights  to  the  boxes  can  alter  the  center  of  gravity  to  simulate  the  effect 
of  equipment  attached  to  the  outside  of  a  helmet.  The  helmet  simulator  has 
been  calibrated  for  weights  of  1.4,  1.8,  2.3,  3.2,  and  4.1  kg.  for  each  of 
the  C.G.  locations  shown  in  Table  2.  Figure  1  illustrates  the  range  of  C.G. 
variations  together  with  definition  of  the  coordinate  axes  by  which  the  C.G. 
locations  are  measured. 

As  shown  in  Figure  5,  a  point  midway  between  the  left  and  right  ear 
canals  has  been  chosen  as  the  origin  of  the  coordinate  axes.  The  helmet 
simulator  is  pictured  in  lateral  (Fig.  6),  oblique  (Fig.  7)  and  frontal 
(Fig.  8)  views.  It  has  been  provided  with  adjustment  to  ensure  that  an 
index  point  on  it  can  be  aligned  with  the  ear  canals,  and  also  with 
independent  adjustment  to  permit  the  suspension  system  to  be  made 
comfortable . 

Eight  headgear  centers-of-gravity  for  three  different  headgear  weights 
(a  total  of  twenty-four  headgear  combinations)  were  evaluated  (as  per  Tables 
3  and  4)  utilizing  the  variable  center-of-gravity  and  variable  weight  helmet 
simulator.  The  "essential  equivalency"  between  the  variable  center-of- 
gravity  and  variable  weight  helmet  simulator  and  selected  headgear  loading 
configurations  has  been  reported  by  Phillips  and  Petrofsky  (1982b). 


EXPERIMENTAL  PROTOCOL 

The  experimental  protocol  may  be  summarized  as  follows: 

Pre-Exercise  MVC 

With  the  subject  seated  in  the  helmet  dynamometer,  the  subject  would 
then  either  perform  a  brief  (3  second)  forward  MVC  (with  EMG  recorded  from 
the  sternocleidomastoid  muscle)  or  a  brief  (3  second)  right  lateral  MVC 
(with  EMG  recorded  simultaneously  from  both  the  posterior  neck  muscles  and 
sternocleidomastoid  neck  muscle).  The  contraction  mode  selected,  would  then 
be  repeated  at  3  minute  intervals  until  3  such  contractions  were  performed. 
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The  strongest  contraction  (highest  strength  and  highest  RMS  amplitude  of  the 
EMG)  would  then  be  taken  as  the  reference  contraction. 

Head  Load inn  Conf iauration  and  Exercise  Duration 

With  the  subject  removed  from  the  isometric  head  dynamometer, 
alternating  right  and  left  lateral  neck  rotations  were  performed  while 
wearing  the  variable  center-of-gravity  and  variable  weight  helmet  simulator 
which  was  set  to  one  of  the  fifteen  headgear  combinations.  The  exercise 
duration  was  30  minutes. 

Post-Exercise  Contractions 

Immediately  upon  completion  of  the  exercise  period,  the  subject 
repositioned  himself  in  the  isometric  head  dynamometer,  and  a  target  tension 
of  70%  of  the  pre-exercise  WVC  was  sustained  (in  the  direction  of  the 
pte-exercise  MVC),  and  held  to  fatigue.  The  duration  of  this  was  called  the 
endurance  time. 

The  order  of  presentation  of  the  direction  of  the  pre-exercise  MVC,  the 
head  loading  configuration,  and  post-exercise  contractions  selected  were  all 
randomized  for  all  of  the  subjects. 

DEFINITION  OF  BOUNDARY  CONDITIONS 

The  limits  of  displacement  along  the  X-  and  Y-axes  are  shown  in  Figure  9 
for  the  helmet  simulator.  The  purpose  of  this  section  is  to  convert  this 
frame  of  reference  into  a  rectangular  system  of  an  equivalent  area.  The  new 
X  and  Y  coordinates  are  the  boundary  conditions  for  the  computer  model. 

The  area  (T)  of  the  shaded  sections  of  Figure  9  is  calculated  as  follows 
T  =  (.5)  (.7)  (1.8)  +  (.5)  (.7)  (1.8) 

+  (2.5)  (.7)  +  (.5)  (.7)  (1.8) 

+  (.5)  (.7)  (1.8) 

T  =  0.63  +  0.63  +  1.75  +  0.63  +  0.63 
T  =  4.27  cm2 

The  total  area  (T1)  is  the  area  of  A  and  B  plus  T: 

T1  ^  (4.3  +  1.8)  (1.8)  +  4.27 
T1  =  10.98  +  4.27  =  15.25  cm2 
Next  define  an  incremental  AX  and  AY  such  that: 

(6.1  +  2AX)  (1.8  +  AY)  =  15.25 
Simplify  by  defining: 

AX  =  AY 

Also,  the  boundaries  are  radial  arcs,  not  straight  lines,  so  that  T  is 
slightly  (12%)  larger  than  previously  calculated: 

T  =  4.76  cm  2 
So  rewriting: 

(6.1+  2AX)  (1.8  +  AX)  =  15.74 
Expanding : 

10.98  +  S.1AX  +  3 . 6 AX  +  2AX  =  15?74 
Rearranging : 

(AX)2  +  4 . 85 (AX)  -  2.38  =  0 


AY  =  AX  =  0.45 


Check: 

(6.1  +  .9)  (1.8  +  .45)  =  15.75 

So  that  the  rectangular  x-y  coordinate  system  is: 

+  X1  =  4.3  +  0.45  =  4.75  cm 

-  X1  =  1.8  +  0.45  =  2.25  cm 
Y1  *  1.8  +0.45  =  2.25  cm 
which  is  shown  in  Figure  10. 

TESTING  OF  ASSUMPTIONS 

Referring  once  again  to  Figure  1,  the  important  physical  parameters  of 
the  system  are:  F,  the  load;  Mjj,  the  moment  with  respect  to  the  X-axis; 
and  My,  the  moment  with  respect  to  the  Y-axis. 

Note  that  at  these  specific  loading  points,  there  is  no  moment  with 
respect  to  the  Z-axis  (i.e.,  they  are  still  parallel  to  it).  That  the 
Z-axis  is  not  considered  to  be  physiologically  significant  is  why  we  chose 
point  CH  for  our  original  study.  In  essence,  the  origin  of  the  three  axes 
("0")  sees  the  same  F,  whether  at  CL  or  at  CH,  and  (of  course)  =  My  =  0. 
This  was  tested  by  observing  whether  forward  and  lateral  contraction 
endurance  times  (for  the  neck  muscles)  with  either  the  CH  or  CL  configuration 
were  similar  (i.e.,  not  statistically  significantly  different).  Furthermore, 
we  tested  whether  the  effects  of  neck  muscle  loading  were  axi-synmetric 
(physiologically).  That  endurance  times  for  C.G.  displacement  along  the  -Y 
axis  (left  side  of  the  head/helmet)  were  similar  to  endurance  times  along 
the  +Y  axis  (right  side  of  the  head/helmet)  is  why  position  LLL  was  selected 
in  our  final  study.  This  was  tested  by  observing  whether  forward  and 
lateral  neck  muscle  contraction  endurance  times  were  similar  for  either  the 
LLL  or  RLL  configuration  (i.e.,  not  statistically  significantly  different). 
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MATHEMATICAL  MBTHODS 

The  model  utilizes  a  three-dimensional  space  defined  by  the  X,  Y  and  F 
axes  (see  Fig.  11).  The  X-axis  defines  the  location  of  the  C.G.  forward 
from  the  system  origin  to  4.75  cm  anterior.  The  Y-axis  defines  the  location 
of  the  C.G.  lateral  from  the  system  origin  to  2.25  cm  right  lateral.  The 
F-axis  defines  the  helmet  configuration  load  from  1.45  kg.  minimum  weight  to 
4.09  kg.  maximum  weight. 

The  boundary  conditions  of  the  three-dimensional  space  are  twelve 
endurance  times  (E^  to  Ej^)*  which  can  be  either  the  forward  contraction  or 
lateral  contraction  endurance  time  to  sustain  an  isometric  neck  muscle 
contraction  at  70%  of  the  muscle's  MVC.  The  boundary  conditions  are 
specified  as  Ey.^.p,  i.e.,  the  endurance  time  after  30  minutes  of  dynamic 
exercise  with  a  helmet  C.G.  located  y  cms  along  the  Y-axis,  x  cms  along  the 
X-axis  and  a  helmet  weight  of  F  kgs. 

The  three-dimensional  system  consists  of  eleven  element  faces  (A  through 
L)  as  shown  at  the  bottom  of  Figure  11-  Each  element  face  requires  a  set  of 
four  subsidiary  equations  (1  to  4)  and  three  primary  equations  (5  to  7)  to 
describe  the  coordinate  (E,  endurance)  point  for  that  face.  Face  A  is  shown 
in  Figure  12. 

Face  A  equations  are: 

Al)  E  =  E  +  (F  -  2.27)  (E  -  E  ) 

AF1  2  1702  3  2 


A2 )  E 

AF2 


E  +  (F  -  2.27)  (E  -  E  ' 
5  T702  6  5 


A3) 


E  =  E  + 
AY1  2 


<  Y)  (E 
TTTS  5 


E  ) 
2 


A4 )  E  +  E  +  (  Y _ )  (E  -  E  ) 

AY2  3  772?  6  3 


A5)  E  =  E  +  (  Y _ )  (E  -  E  ) 

AF  AF1  77775  AF2  AF1 


A6)  E 

AY 


E  +  (F  -  2.27)  (E 
AY1  1702  AY2 


K  ) 
AY1 


A7)  E  =  E  +  E 

A  AF  AY 

- 2 - 

Before  we  proceed  to  FACES  B,  C,  D,  F,  G,  H,  I,  J,  K ,  L  (in  a  similn-  manner), 
we  repeat  the  process  for  the  standard  deviations  (SD’s). 
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The  boundary  conditions  for  the  three  dimensional  space  can  also  be  the 
twelve  standard  deviations  (SD^  to  SDj^)  of  the  mean  endurance  times 
described  previously.  This  is  shown  in  Fig.  13  as  a  system  of  SD’s 
(S ■ D. y i x • F> ■  T*'e  same  eleven  face  elements  (AA  to  LL)  are  present  as  with 

the  endurance  times.  Each  element  face  is  described  by  the  same  set  of 
seven  equations.  Face  AA  is  shown  in  Fig.  14. 

Face  AA  equations  are: 


AA1) 

to 

o 

11 

SD  +■ 

(F  -  2. 

27) 

(SD  -  SD  ) 

AF1 

2 

1.57 

3  2 

AA2) 

SD 

=  SD  + 

(F  -  2. 

27) 

(SD  -  SD  ) 

AF2 

5 

1.87 

6  5 

AA3) 

SD 

=  SD  + 

(  Y  ) 

(SD 

-  3D  ) 

AF2 

2 

7775 

5 

2 

AA4 ) 

SD 

=  SD  +- 

(  Y  ) 

(SD 

-  SD  ) 

AY2 

3 

7775 

6 

3 

AA5) 

SD  ■ 

SD  + 

(Y  ) 

(SD 

-  SD  ) 

AF 

AF1 

2^5 

AF2  AF1 

AA6 ) 

SD  = 

SD  + 

(F  -  2 

.27) 

(SD  -  SD  ) 

AY 

AY1 

1.82 

AY2  AY1 

AA7) 

SD  = 

SD  + 

SD 

A 

AF 

AY 

- 2 

We  now  proceed  to  FACES  B  through  L,  but  now  computing  an  SD  for  each 
face  (i.e.,  replace  E  with  SD) .  Results  are  given  in  Appendix  A. 

For  any  load  (F)  between  2.27  kg.  and  4.09  kg.,  a  particular  endurance 
time  E  can  be  interpolated  using  equations  Ml  to  M4  as  shown  in  Figure  15. 

For  any  load  (K)  between  1.45  kgs.  and  2.26  kgs.,  a  particular  endurance 
time  E  can  be  interpolated  using  equations  M5  to  M8  as  shown  ir.  Figure  16. 

Also  recall  since  we  have  computed  the  SD's  for  each  of  the  11  faces: 


MH1)  SD  =  SD  +  (  X  )  (SD  -  SD  ) 
AF  A  Z~J5  F  A 


MM2)  SD  =  SD  +  ( _ Y  )  (SD  -  SD  ) 

CH  C  7725  H  C 
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MM3) 

SD 

=  SD  x  (F  -  2. 

27)  (SD  -  SD  ) 

KJ 

X  T7TT2 

J  K 

MM4) 

SD  = 

SD  +  SD  +  SD 

AF  CH  KJ 

- 3 - 

MM5) 

SD 

=  SD  +  (  X  ) 

(SD  -  SD  ) 

BG 

B  *775 

G  B 

MM6) 

SD 

a  SD  +  (  Y  ) 

(SD  -  SD  ) 

DI 

D  7775 

I  D 

MM7) 

SD 

*  SD  +  (F  -  1. 

45)  (SD  -  SD  ) 

LK 

L  0752 

K  L 

MM  8  / 

SD  = 

SD  +  SD  +  SD 

BG  DI  LK 
- 3 - 


The  other  half  of  the  model  utilizes  a  three-dimensional  space  defined 
by  the  -X,  Y,  and  F  axes  (see  Figure  17).  The  -X  axis  (or  the  W-axis) 
defines  the  location  of  the  C.G.  backward  from  the  system  origin  to  2.25  cm 
rearward.  The  Y-axis  and  F-axis  are  as  previously  defined.  The  boundary 
conditions  of  the  three-dimensional  space  are  the  four  previous  endurance 
tiraea  (Elt  E3,  E^  and  E6)  and  eight  more  endurance  times  (E13  to  E18). 

The  three-dimensional  system  consists  of  eleven  element  faces  (A  end  B 
previously,  as  well  as  N  to  V)  as  shown  at  the  bottom  of  Figure  17.  Each 
element  face  requires  the  same  seven  equations  to  describe  the  coordinate 
(E,  endurance)  point  for  that  face.  Note  that  Faces  A  and  B  were  previously 
defined.  Equations  for  Faces  N  to  V  are  given  in  Appendix  B.  A  system  of 
SB* 6  (SDy.y.p)  can  also  be  defined  for  face  elements  (NN  to  W). 

For  any  load  (F)  between  2.27  kg.  and  4.09  kg.,  a  particular  endurance 
time  E  can  be  interpolated  vjing  equations  W1  to  W4  as  shown  in  Figure  18. 

For  any  load  (F)  between  1.45  kg.  and  2.26  kg.,  a  particular  endurance 
time  E  can  be  interpolated  using  equations  W5  to  W8  as  shown  in  Fig.  19. 

Also  recall  since  we  have  computed  the  SD's  for  each  of  the  11  faces: 


WW1)  SD  =  SD  +  (  W _ )  (SD  -  SD  ) 

AP  A  7775  P  A 


WW2)  SD  =  SD  +  (  Y _ )  (SD  -  SD  ) 

RM  F  7775  N  R 


WW3)  SD  =  SD  +  (F  -  2.27)  (SD  -  SD  ) 
UT  U  1752  T  U 
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WW4)  SO  a  SO  -I-  so  +  SO 


AP  RN 

UT 

3 

WW5) 

SD 

SD 

+ 

<  W  ) 

(SD  -  SD  ) 

BQ 

B 

272F 

Q  B 

WW6) 

SD 

SD 

-*■ 

(  Y  ) 

(SD  -  SD  ) 

SO 

S 

772T 

0  S 

WW7) 

SD 

SD 

+ 

(F  -  1. 

45)  (SD  -  SD  ) 

VU 

V 

U  V 

WW8)  SD  =  SD  +  SO  +  SD 

BO  SO  VU 

— 54 - !r - 


We  now  modify  the  model  to  account  for  specific  endurance  timer 
associated  with  forward  or  lateral  contractions.  Suppose  for  through  Bjq, 
these  are  the  endurance  times  for  forward  contractions.  Then  define: 

B  =  EF 

and  substitute  this  term  in  all  the  preceding  equations,  i.e.,  A1  through  W8. 

Suppose  now  that  E^  through  E^g  represent  the  respective  endurance 
times  for  lateral  contractions: 

Then  define: 


E  =  EL 

and  substitute  in  all  equations  (A1  through  W8) . 

Consequently,  all  equations  (A1  through  W8)  must  be  written  twice: 
For  example,  regarding  forward  contractions: 


Al)  EF  =  BP  +  (F  -  2.27)  (EF  -  KF  >  to  .  . 
AF1  2  17^2  3  2 


W8)  EF 


EF 


+  EF 

SO  + 

- j - 


EF 

_ VU 


For  example,  regarding  lateral  contractions: 

Al)  EL  =  EL  +  (F  -  2.27)  (EL  +  EL  )  to  .  .  . 
AF1  2  Y7W2  3  2 

W8)  EL  =  EL  +  EL  +  EL 

BQ  SO  _ VU 
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The  process  can  be  easily  repeated  for  the  respective  equations  that 
solve  for  the  standard  deviation  (SO). 

For  forward  contractions,  define: 

SD  *  SF 

and  substitute  this  term  in  all  the  parallel  equations  (AA1  through  WW8) . 
For  lateral  contractions,  define: 

SD  =  SL 

and  substitute  this  term  in  all  the  parallel  equations  (AA1  through  WW8) . 

The  computer  program  must  be  directed  to  solve  for  E  and  SD  in  one  of 
four  quadrants  as  shown  in  Figure  20. 

This  can  be  done  by  BASIC  statements  "IF. . .THEN. . . ”  which  process  the 
input  data  and  direct  the  computer  to  the  correct  quadrant  (i.e. ,  sets  of 
equations) . 

IMP  ITT  DATA  (Range  of  Values): 

F  (Load):  1.45  kg  to  4.09  kg 

X  (X-Axis):  -2.25  cm  to  4.75  cm 

YIN  (Y-Axis):  -2.25  cm  to  2.25  cm 

DATA  PREPROCESSING: 


60  Y  =  ABS  (YIN) 


65 

IF 

F 

< 

1.45 

THEN 

1000 

70 

IF 

F 

> 

4.09 

THBN 

1000 

75 

IF 

Y 

> 

2.25 

THEN 

1000 

80 

IF 

X 

> 

4.75 

THEN 

1000 

85 

IF 

X 

< 

-2.25 

THEN 

1000 

1000  PRINT  "DATA  OUT  OF  BOUNDS" 
1005  END 


DATA  ROUTING 

This  is  a  series  of  tests  to  partition  the  data  in  the  appropriate 
quadrant . 

100  IF  F  <  2.27  THEN  300 

105  IF  X  <  0  THEN  500 

110  REM:  Equations  for  Quadrant  A: 

A1  -  A7 ,  Cl  -  C7 ,  FI  -  F7 ,  HI  -  H7, 

J1  -  J7,  K1  -  K7,  Ml  -  M4; 

184  AA1  -  AA7,  CC1  -  CC7,  FF1  -  FF7, 

Equations  HH1  -  HH7 ,  JJL  -  JJ7 ,  KK1  -  KK7, 

MM1  -  MM4 ;  (Equations  written  for  EF,  EL,  SF,  SL) . 

290  GOTO  900 

300  IF  X  <  0  THEN  700 

305  REM:  Equations  for  Quadrant  B: 

B1  -  B7 ,  D1  -  D7 ,  Cl  -  G7 ,  II  -  17, 

184  K1  -  K7,  LI  -  L7,  M5  -  M8; 
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Equations 

BB1  -  BB7, 

DD1  -  DD7 ,  GG1  -  GG7 , 

III  -  117 , 

KK1  -  KK7,  LL1  -  LL7 , 

MM5  -  MM8; 

(Equations  written  for  EF,  EL,  SF,  SL) 

490 

GOTO 

900 

500 

U  a  ■ 

-X 

505 

REM: 

Equations 

for  Quadrant  C: 

A1  -  A7 ,  Ml 

-  N7,  PI  -  P7,  R1  -  R7 , 

184  T1  -  T7,  U1  -  U7,  W1  -  W4 ; 

Equations  AA1  -  AA7 ,  NN1  -  NN7 ,  PP1  -  PP7 , 

HR1  -  RR7 ,  TTJ  -  TT7,  UU1  -  UU7, 

VWl  -  WW4;  (Equations  written  for  EF,  EL,  SF,  SL) 

690  GOTO  fOO 

700  W  =  -X 

705  RBM:  Equations  for  Quadrant  D: 

B1  -  B7,  01-07,  Q1  -  Q7 ,  SI  -  S7 , 

184  U1  -  U7,  VI  -  V7,  W5  -  W8; 

Equations  BB1  -  BB7 ,  001-007,  QQ1  -  QQ7 , 

SSI  -  SS7 ,  UU1  -  UU7 ,  VV1  -  W7 , 

WW5  -  WW8;  (Equations  written  for  BF,  EL,  SF,  SL) 

900  PRINT  ‘’RESULTS  ARE .  .  . " 

995  GO  TO  1005 

The  overall  computer  flow  chart  which  results  is  shown  in  Figure  21. 

RESULTS 

The  results  are  summarized  into  four  categories:  the  endurance  times 
(and  their  standard  deviations)  which  represent  the  boundary  conditions,  the 
vertical  endurance  time,  the  horizontal  endurance  time,  and  the  final 
computer  model. 

ENDURANCE  TIME  (BOUNDARY  CONDITIONS) 

The  mean  endurance  times  (and  standard  deviations)  which  were 
experimentally  determined  and  used  as  the  boundary  conditions  for  the 
mathematical  model  are  presented  in  Table  5.  The  stored  constants  used  in 
the  computer  program  (top  of  Figure  21),  EFj  to  KF^g  and  SFj  to  SFjg,  are 
those  values  presented  in  Table  5A.  The  other  stored  constants,  ELj  to  BL18 
and  SLjto  SLjg,  are  those  values  shown  in  Table  5B.  This  is  summarized  in 
Table  6 . 

ENDURANCE  TIMB  (VERTICAL) 

The  endurance  times  (and  standard  deviations)  for  helmet  weights  of  1.45 
kg.,  2.27  kg.  and  4.09  kg.  when  the  C.G.  was  center-high  (CH)  compared  to 
center-low  (CL)  are  presented  for  forward  isometric  neck  muscle  contractions 
in  Table  7A. 

The  endurance  times  (and  standard  deviations)  for  the  same  three  helmet 
weights  and  same  two  C.G.’s  are  presented  for  lateral  isometric  neck  muscle 
contractions  in  Table  7B. 

For  both  tables,  paired  T-tests  between  the  two  C.G.  locations  (for  each 
of  the  three  helmet  weights)  showed  no  statistically  significant  difference. 
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ENDURANCE  TIKE  (LATERAL) 

The  endurance  times  (and  standard  deviations)  for  helmet  weights  of  1.45 
kg.,  2.27  kg.  and  4.07  kg.  when  the  C.G.  was  left-lateral-low  (LLL)  compared 
to  right-lateral-low  (RLL)  are  presented  for  forward  isometric  neck  muscle 
contractions  in  Table  8A. 

The  endurance  times  (and  standard  deviations)  for  the  same  thret.  helmet 
weights  and  same  two  C.G.'s  are  presented  for  lateral  isometric  neck  muscle 
contractions  in  Table  8B. 

For  both  tables,  paired  T-tests  between  the  two  C.G.  locations  were  not 
significantly  different  at  1.45  kg.  and  2.27  kg.  of  helmet  weight.  The  two 
C.G.  locations  were  significantly  different  (p  <  -05)  at  the  4.09  kg. 
helmet  weight. 

COMPUTER  MODEL 

The  final  computer  model  was  constructed  utilizing  the  stored  constants 
EPi  -  EF10,  SFx  -  SF18,  BLj  -  EF18  and  SLx  -  SL18  (as  per  Table  6).  EFj  - 
EFis  were  substituted  into  Eqs.  A1  to  W8  and  ELj  -  EL^g  were  also  substituted 
in  Eqs.  A1  to  W8  as  per  the  Mathematical  Methods  (see  METHODS  AND  MATERIALS 
section).  SFi  -  SFj8  and  SL^  -  SL^g  were  substituted  into  Eqs.  AA1  to  WW8 . 
The  program  is  user  interactive  and  a  complete  listing  is  given  in  Appendix 

m. 

The  user  enters  the  helmet  weight  (in  kilograms),  the  X-axis  coordinate 
of  the  C.G.  and  the  Y-axis  coordinate  of  the  C.G.  (in  centimeters). 

The  program  will  then  print  out: 

(a)  the  mean  endurance  time  of  a  forward  isometric  neck  muscle 
contraction  (sustained  at  70%  of  the  MVC)  in  seconds; 

(b)  the  standard  deviation  of  (a); 

(c)  the  mean  endurance  time  of  a  lateral  isometric  neck  muscle 
contraction  (sustained  at  70%  of  the  MVC)  in  seconds; 

(d)  the  standard  deviation  of  (c). 

A  special  significance  test  program  has  also  been  written  and  is  listed 
in  Appendix  IV.  If  a  second  helmet  weight  and  C.G.  combination  is  also 
evaluated,  the  significance  test  program  will  test  as  to  whether  the  first 
helmet  weight  and  C.G  loading  combination  is  significantly  diferent  from  the 
second  helmet  weight  and  C.G.  loading  combination  with  respect  to  neck 
muscle  endurance  times. 


DISCUSSION  AND  CONCLUSIONS 

A  computer  model  of  neck  muscle  endurance  and  fatigue  as  a  function  of 
helmet  loading  has  been  developed.  The  final  model  consists  of  over  700 
equations  and  has  been  formatted  to  run  on  an  Apple  11+  and/or  Apple  lie 
microcomputer  with  at  least  48K  of  memory. 

Our  objective  has  been  to  make  the  model  "definitive,"  but  the  final 
model  deviates  from  this  desired  objective  in  several  respects.  An 
appreciation  of  those  limitations  will  allow  the  user  to  more  reasonably 
interpret  results  he/she  may  obtain  from  the  model. 
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First,  the  model  it  valid  only  over  a  limited  range  of  weight  and  C.G. 
locations.  Specifically,  helmet  weight  must  be  between  1.45  kg  (3.2  lbs) 
and  4.09  kg  (9.0  lbs),  the  X-axis  displacement  must  be  between  -2.25  cm  and 
4.75  cm  and  the  Y-axis  displacement  must  be  between  -2.25  cm  and  +2.25  cm. 

Any  conditions  outside  these  limits  will  cause  the  program  to  print  a  "DATA 
OUT  OF  BOUNDS"  statement. 

Second,  the  model  has  been  defined  for  a  three-dimensional  space  based 
upon  defined  boundary  conditions  (see  Methods  and  Materials.  Mathematical 
Methods) .  It  assumes  that  any  helmet  weight  and  C.G.  location  that  falls 
within  the  three-dimensional  space  can  be  linearly  interpolated  between  the 
various  boundary  conditions.  In  reality,  the  model  is  "piece-wise  linear" 
between  helmet  weights  of  1.45  kg  to  2.26  kg  2.27  kg  to  4.09  kg.  It  is 

also  "piece-wise"  linear  for  X-axis  displacements  between  -2.25  cm  to  0  era 
and  O  cm  to  4.75  cm.  Finally,  it  is  "piece-wise"  linear  for  Y-axis 
displacements  between  -2.25  cm  to  0  cm  and  0  cm  to  +2.25  cm. 

Third,  the  assumption  of  insensitivity  to  vertical  loading  has  been 
established,  but  only  for  the  0,0  X-Y  coordinate.  This  assumption  was  not 
tested  at  any  of  the  other  boundary  conditions.  Furthermore,  flight 
situations  in  which  all  G  forces  are  not  vertical,  (i.e.,  where  there  is 
some  lateral  loading  such  as  Gj  or  Gy  forces)  may  make  the  vertical 
insensitivity  assumption  invalid. 

Fourth,  the  assumption  of  axi-symmetry  (i.e.,  data  aquired  for  right 
sided  head  loading  is  also  valid  for  left  sided  head  loading)  has  only  been 
partially  validated.  The  assumption  appears  valid  for  light-to-moderate 
helmet  weights  (i.e.,  1.45  kg  to  2.27  kg)  but  again  was  only  tested  at  the 
0,0  X-Y  coordinate  system.  This  assumption  was  not  tested  at  any  other 
boundary  condition.  Furthermore,  for  heavy  helmet  weights  (4.09  kg)  the 
assumption  was  not  validated.  Results  predicted  by  the  model  for  heavier 
weight  helmets  with  left  sided  C.G.  locations  must  be  viewed  with  discretion. 

Fifth  and  finally,  laboratory  based  studies  which  predict  isometric  neck 
muscle  endurance  times  (either  forward  or  lateral)  as  a  function  of  headgear 
loading  cannot  be  extrapolated  to  predict  "operational  mission  endurance 
times"  (OMET).  OMET's  are  a  complex  function  of  operator  workload,  heat, 
noise,  vibration,  etc. 

Therefore,  OMET's  are  multi-factorial  problems  which  do  not  easily  lend 
themselves  to  mathematical  analysis.  It  was  the  inability  of  our  group  to 
translate  our  laboratory  based  isometric  endurance  times  into  meaningful 
OHET’s  that  was  a  major  factor  in  not  pursuing  this  project  further.  We 
remain  open  to  suggestions  and  would  be  pleased  to  submit  a  proposal  in  this 
area  if  warranted. 


Subiect 

Age  (yrs . ) 

Height  (cm) 

1 

22 

130 

2 

28 

175 

3 

22 

193 

4 

20 

188 

5 

20 

185 

6 

19 

168 

Weight  (kg.)  Meek  Cir.  (cm) 


77 

37 

85 

37 

79 

38 

73 

37 

70 

36 

61 

36 

TABLE  1.  GENERAL  CHARACTERISTICS  OF  THE  SUBJECTS 
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C.G.  Location 
Description* 


Central 


Maximum  Height 

0 

0 

8.0 

Medium  Height 

0 

0 

4.0 

Low  Height 

0 

0 

0 

Forward 

Maximum  Height 

5.0 

0 

8.0 

Medium  Height 

S.O 

0 

4.0 

Low  Height 

5.0 

0 

0 

Af  tward 

Maximum  Height 

-2.5 

0 

8.0 

Medium  Height 

-2.5 

0 

4.0 

Low  Height 

-2.5 

0 

0 

Central 

Left,  Maximum  Height 

0 

2.5 

8.0 

Left,  Medium  Height 

0 

2.5 

4.0 

Left,  Low  Height 

0 

2.5 

0 

Central 

Right,  Maximum  HeighL 

0 

-2.5 

8.0 

Right,  Medium  Height 

0 

-2.5 

4.0 

Right,  Low  Height 

0 

-2.5 

0 

Forward 

Left,  Maximum  Height 

4.3 

1.8 

8.0 

Left,  Medium  Height 

4.3 

1.8 

4.0 

Left,  Low  Height 

4.3 

1.8 

0 

Forward 

Right,  Maximum  Height 

4.3 

-1.8 

8.0 

Fight,  Medium  Height 

4.3 

-1.8 

4.0 

Right,  Low  Height 

4.3 

-1.8 

0 

*  C.G.  locations  for  total  weights  of  less  than  1.8  kg.  may  differ  from 
values  shown  in  this  table. 

**  Displacement  from  head  and  neck  c.g.  (axis  directions  as  defined  in 
Figure  1). 

TABLE  2.  SPECIFIC  C.G.  LOCATIONS  FOR  WHICH  THE  HELMET 
SIMULATOR  IS  CALIBRATED 
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Til 


/WVW.s 


C.G. -Weight 
Description 

Weight 

(kgs) 

X 

Displacement 

Y 

(cm)* 

2 

CL 

Center-Low-3 

1.45 

0 

0 

0 

Center-Low-5 

2.27 

0 

0 

0 

Centor-Low-9 

4.09 

0 

0 

0 

CH 

Center-High-3 

1.45 

0 

0 

8.0 

Center-High-5 

2.27 

0 

0 

8.0 

Center-High-9 

4.09 

0 

0 

8.0 

FL 

Forward-Low-3 

1.45 

5.0 

0 

0 

Forward-Low-5 

2.27 

5.0 

0 

0 

Forward-Low-9 

4.09 

5.0 

0 

0 

LRL 

Lat-Right-Low- 3 

1.45 

0 

-2.5 

0 

Lat-Right-Low-5 

2.27 

0 

-2.5 

0 

Lat-Right-Low-9 

4.09 

0 

-2.5 

0 

AL 

Af  terward-Low-3 

1.45 

-2.5 

0 

0 

Afterward-Low- 5 

2.27 

-2.5 

0 

0 

Af t e rwa rd - Low- 9 

4.09 

-2.5 

0 

0 

*  Displacement  from  the  head/neck  center-of-gravity :  axis  directions  as  per 
Fig.  1,  and  displacement  distances  as  per  lable  2. 


TABLE  3.  FIFTEEN  HEADGEAR  COMBINATIONS  ORIGINALLY  EVALUATED 
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C . G . -Weight 
Description 

Weight 

(Kgs) 

X 

Displacement 

Y 

(cm)* 

Z 

FBL 

For-Right~Low-3 

1.45 

4.3 

-1.8 

0 

For-Right-Low-5 

2.27 

4.3 

-1.8 

0 

For-Right-Low-9 

4.09 

4.3 

-1.8 

0 

ARL 

Af t-Rlght-Low-3 

1.45 

-1.8 

-1.8 

0 

Aft-Right-Low-5 

2.27 

-1.8 

-1.8 

0 

Aft-Right-Low- 9 

4.09 

-1.8 

-1.8 

0 

LLL 

Lat-Lef t-Low-3 

1.45 

0 

2.5 

0 

Lat-Lef t-Low-5 

2.27 

0 

2.5 

0 

Lat-Lef t- Low-9 

4.09 

0 

2.5 

0 

Same  as  footnote  (Table  3.) 


TABLE  A.  NINE  ADDITIONAL  HEADGEAR  COMBINATIONS  EVALUATED 
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A.  Forward  Contractions 


Wt. 


*C .  G . 


1.45  kg 


2.27  kg. 


4.09  kg. 


FRL 

49 

+ 

19 

37 

+ 

10 

41 

+ 

9 

FL 

86 

+ 

49 

65 

+ 

26 

56 

+ 

24 

LRL 

75 

jr 

57 

66 

+ 

21 

96 

♦ 

61 

CL 

98 

+ 

34 

74 

+ 

43 

68 

+ 

34 

ARL 

43 

+ 

12 

46 

+ 

16 

45 

+ 

13 

AL 

86 

+ 

34 

64 

+ 

19 

83 

+ 

20 

B.  Lateral  Contractions 


Wt. 


*C.G. 


1.45  kg. 


2.27  kg. 


4.09  kg. 


FRL 

105 

t 

25 

102 

+ 

26 

122 

+ 

49 

FL 

148 

+ 

53 

89 

+ 

42 

62 

+ 

35 

LRL 

156 

± 

59 

81 

+ 

40 

72 

+ 

40 

CL 

106 

+ 

88 

83 

+ 

33 

86 

+ 

38 

ARL 

82 

+ 

31 

79 

47 

94 

+ 

31 

AL 

91 

+ 

38 

73 

47 

129 

+ 

69 

See  "CG-Weight  Description”  Column  (Table  3). 


TABLE  5.  ENDURANCE  TIME:  BOUNDARY  CONDITIONS  (Secs) 
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.  t  -  ■  a  TV  .  ft  _  Hi.  Aaf.aP.a 


n  a*  r.  r 


*_A 


CONSTANT 

VALUE 

CONSTANT 

VALUE 

CONSTANT 

VALUE 

CONSTANT 

VALUE 

BFi 

98 

SF1 

34 

EL1 

106 

SL1 

88 

bf2 

74 

sf2 

43 

EL2 

83 

SL2 

33 

BF3 

68 

sf3 

34 

el3 

86 

SL3 

38 

BF4 

75 

SF4 

57 

BL4 

156 

SL4 

59 

bf5 

66 

sf5 

21 

el5 

81 

sl5 

40 

BF6 

96 

SF6 

61 

EL6 

72 

SL6 

40 

BF? 

49 

sf? 

19 

EL7 

105 

sl? 

25 

EF8 

37 

SF8 

10 

EL8 

102 

SL8 

26 

KF9 

41 

SF9 

9 

EL9 

122 

SL9 

49 

BFio 

86 

srio 

49 

elio 

148 

slio 

53 

BF11 

65 

SF11 

26 

EL11 

89 

SL11 

42 

SFi2 

56 

SF,, 

12 

24 

BL 

12 

62 

SL 

12 

35 

KF13 

43 

SF 

13 

12 

EL 

13 

82 

SL10 

13 

31 

BF14 

46 

SF 

14 

16 

EL 

14 

79 

SL 

14 

47 

BF 

15 

45 

SF 

15 

13 

EL 

15 

94 

SL15 

31 

EF 

16 

86 

SF. 

16 

34 

EL 

16 

91 

SL16 

38 

BF 

17 

64 

SF17 

19 

EL17 

73 

SL 

17 

47 

EF 

18 

83 

SF18 

20 

E4s 

129 

SL 

18 

69 

TABLE  6.  IDENTIFICATION  OF  STORED  CONSTANTS 
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A.  Forward  Contractions 


*C.G. 

1.45  kg. 

2.27  kg. 

4.09  kg. 

CH 

71  +  49 

82  +  59 

71  +  55 

CL 

98  +  34 

74  +  43 

68  +  34 

p-level 

N.S. 

N.S. 

N.S. 

B.  Lateral 

Contractions 

Wt. 

*C.G. 

1.45  kg. 

2.27  kg. 

4.09  kg. 

CH 

82  +  56 

72  +  35 

61  +  39 

CL 

106  +  88 

83  ±  33 

86  +  38 

p-level 

N.S. 

N.S. 

N.S. 

*  See  footnote  (Table  5). 


TABLE  7.  EMDURANCE  TIME:  VERTICAL  (Secs) 
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A.  Forward  Contractions 


Wt. 


*C.G . 

1.45  kg- 

2.27  kg. 

4.09  kg. 

LLL 

46  ±  18 

47  +  19 

43  +  10 

RLL 

75  +  57 

66  ±  21 

96  +  61 

p-level 

N.S. 

N.S. 

p  <  .05 

B.  Lateral 

Contractions 

Wt. 

*C.G. 

1.45  kg. 

2.27  kg. 

4.09  kg. 

LLL 

108  +  50 

110  +  27 

116  +  35 

RLL 

156  +  59 

81  +  40 

72  +  40 

p-level 

N.S. 

N.S. 

p  <  .05 

*  See  footnote  (Table  7). 


TABLE  8.  ENDURANCE  TIME:  LATERAL  (Secs) 
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FIGURE  4.  SUBJECT  USING  THE  HEAD  DYNAMOMETER:  LATERAL  VIEW 
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FIGURE  6.  HELMET  SIMULATOR:  LATERAL  VIEW 
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FIGURE  7.  HELMET 


FIGURE  9 


New  Coordinate  System 
(See  Fig  10) 


T=  Cross- Hatched  Area 


SCHEMATIC  OF  THE  ACTUAL  HELMET  SIMULATOR  COORDINATE  SYSTEM 
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Three  Dimensional  System:  (Ey,x,f) 


~t<E  0,4.75.4.1 

'  \<E12) 
I  V. 


Simplified:  (see  also  SD’s) 


11.  THE  THREE  DIMENSIONAL  SPACE  SYSTEM  (E 


Face  A: 


EAF1  ♦ — 
I 


♦ 


FIGURE  12.  FACE  "A"  ELEMENT 
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Three  Dimensional  System  (SD’s):  SDy.x.f 


(SO  4)  (SD7) 


PIGURE  13.  THE  THREE  DIMENSIONAL  SPACE  SYSTEM  (SD 

7  »xt  1 
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44 


For  F«  1.45  to  2.26  Kg 


M5) 

m 

CD 

O 

i  *  i 

-  Eg  ♦  \4.75y 

1  <eg-eb> 

M6) 

edi 

(  y  1 
•  Ed  +  \2.25) 

1  <ErED> 

M7) 

elk 

/  F-1.45\ 

•  El*V  0.82  /(Ek-El) 

M8) 

E  = 

eBG  +  eDI  +  eLK 

3 

FIGURE  16.  QUADRANT  ’B'  SPACE  SYSTEM  AMD  EQUATIONS 


4b 


Three  Dimensional  System:  (Ey.w.f) 


FIGURE  17.  THE  THREE  DIMENSIONAL  SPACE  SYSTEM  (E 


For  F=1 .45  to  2.26  Kg 


W5) 

W6) 

W7) 

W8) 


EVU  "  EV  +  '  082  '  (Eu'  EV) 

P  «  eBQ  +  Eso  +  Eyu 
3 


FIGURE  19.  QUADRANT  ’D'  SPACE  SYSTEM  AND  EQUATIONS 
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FIGURE  20.  SUMMARY  OF  THE  FOUR  QUADRANT  SYSTBH 


Quadrant  B 


9 


FIGURE  21.  FLOW  DIAGRAM  OF  THB  COMPUTER  PROGRAM 
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1812  I  N°UT  DO  YOU  WANT  A  HARDCOPY7  ENTER  Y  (YES)  OP  N(NC>  "|H» 
1817.  IF  A*  =»  " N "  GOTO  4005 

1814  IF  At  =  "/"  THEN  F-R»  1:  G0SU6  2070 

1815  F  Ptt  o 
1616  GOTO  4605 

1270  T.LM  FT  306  7  3  F  GR  GUAOA'*Fti  A 

1575  HOMe 

184U  POINT  ‘  PL DUi  15  FOR  THt  FOLLOWING  INPUTS  .  P‘F.1  P’RINT 

1845  FRINT  “LOAD  >«  ;F>  F'RINT 

1650  PRINT  “X-AXIS  •  "|Xl  PRINT 

1855  PRINT  “Y-AXIS  -  “|Vli  PRINT 

1659  PRINT  “  **  RESULTS  •  *"«  PRINT 

184,0  PRINT  ENDUP  ANCL  F  OF<wARD  -  "  |  M4 :  PRINT 

1665  PRINT  ‘  ST  0.  DEVI  AT  ION  FORWARD  *  "  1  Ml. :  FRt.NI 

507u  FRIN1  "ENPUTyANCE  LATERAL  «  '  |  MD  i  PRINT 

1975  FhlNl  '  STD.  DEVIF.T  ION  LATERAL.  -  ‘jMTl  FRINT 

1880  Rt  tl,pn 

16V0  RF.M  RESULTS  FOR  CUADRANT  h 

1 895  HOME 

1YUO  HR  IN!  THE  RESULTS,  FOR  THF  FOLLOWING  INPUTS  AR'.  i  I  R I  NT 
1905  PF'INT  "LOAD  »  "  |  F  i  PRINT 
:9|T.  FPInT  "X-AXIS,  -  ",Xi  PRINT 
1915  PRINT  "Y-AXIS  •  ‘ | Y I l  PRINT 


t.-l  3 


1920 

1923 

1930 

1933 

1940 

1943 

1935 

1960 

1965 

1970 

1975 

1980 

1903 

1990 

1993 

2000 

2003 

2010 

2020 

2023 

2030 

2033 

2040 

2043 

2050 

2053 

2060 

2063 

2070 

2073 

4000 

4003 

4006 

4010 

4013 

4020 

9999 


PR  I  NT  "  ••  RESULTS  •  «"|  PRINT 

PRINT  "ENDURANCE  EORW^O  -  - 1  <J»i  PRINT 
PRINT  '"STD.  DEVIATION  FORWARD  -  “|HP»  PRINT 
PRINT  “ENDURANCE  LATERAL  •  "|HHi  PRINT 
PRINT  "STD.  DEVIATION  LATERAL  -  “|MX|  PRINT 
RETURN 

REM  RESULTS  FOR  QUADRANT  C 
HOME 

PRINT  “THE  RESULTS  FOR  THE  FOLLOW I NG  INPUTS  ARE" t  PRINT 

PRINT  "LOAD  -  “Fi  PRINT 

PRINT  "X-AXIS  -  "|Xi  PRINT 

PRINT  "Y-AXIB  “  “ | Y  f  t  PRINT 

PRINT  "  *«  RESULTS  **"i  PRINT 

PRINT  "ENDURANCE  FORWARD  ■  *)W4i  PRINT 

PRINT  "STD. DEVIATION  FORWARD  -  "|WLi  PRINT 

PRINT  “ENDURANCE  LATERAL  -  "|WD|  PRINT 

PRINT  "STD. DEVIATION  LATERAL  -  "|WTi  PRINT 

RETURN 

REM  RESULTS  FOR  QUADRANT  D 
HOME 

PRINT  "  THE  RESULTS  FOR  THE  FOLLOWING  INPUTS  ARE'"i  PRINT 

PRINT  “LOAD  »  “ | F I  PRINT 

PRINT  "X-AXIS  -  "|Xl  PRINT 

PRINT  “Y-AXIS  «  “ | YT t  PRINT 

PRINT  "  ••  RESULTS  ••"'i  PRINT 

PRINT  “ENDURANCE  FORWARD  ■  ""  J  W8l  PRINT 

PRINT  "STD. DEVIAT ION  FORWARD  -  “tWPi  PR  I N r 

PRINT  "ENDURANCE  LATERAL  «  " | WHl  PRINT 

PRINT  "STD. DEVIATION  LATERAL  -  "jWX«  PRINT 

RETURN 

PRINT  "DATA  OUT  OF  BOUNDS  "  I  PRINT 
PRINT  l  PRINT  "TO  REENTER  DATA" 

PRINT  "OR  TO  ENTER  NEW  DATA  TYPE  1  "I  PRINT 

INPUT  TO  END  PROGRAM  TYPE  0  "' |  2 

IF  I  =  1  GOTO  13V 

IF  2  a  0  GOTO  9999 

END 
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APPENDIX  D 


COMPUTER  LISTING  (TBB  TBST  PROGRAM) 


3  MEN  t  AfWV  WELCH  PROGRAM 

10  PRINT  “ONE-TAILED  WELCH  TEST  PROGRAM" 

20  PRINT  “NOTE  t  THIS  PROGRAM  TESTS  FOR  WHETHER  A  LARGER  (SMALLER)  ENDURAN 
CE  TIME  FOR  ONE  HEADGEAR  CONFIGURATION" 

25  PRINT  "IS  SIGNIFICANTLY  LARGER  (SMALLER)  COMPARED  TO  A  SECOND  HEADGEAR 
CONFIGURATION.  ■' 

30  REM  «  C.A.P.  27— MAY-8*1 . 

50  PRINT 

105  INPUT  "WHAT  IS  THE  MEAN  VALUE  OF  THE  FORWARD  ENDURANCE  TIME  FOR  THE  F 
IR9T  HEADGEAR  CONFIGURATION  "sAl 
107  PRINT 

llO  iNPUT  “WHAT  IS  THE  STANDARD  DEVIATION  OF  THE  FORWARD  ENDURANCE  TIME  F 
OR  THE  FIRST  HEADGEAR  CONFIGURATION  1 01 
115  VI  «  Q1  •  01 
122  PRINT 

125  INPUT  “WHAT  IS  THE  MEAN  VALUE  OF  THE  LATERAL  ENDURANCE  TIME  FOR  THE  F 
IRST  HEADGEAR  CONFIGURATION  ?"SA3 
128  PRINT 

130  INPUT  “WHAT  IS  THE  STANDARD  DEVIATION  OF  THE  LATERAL  ENDURANCE  TIME  F 
OR  THE  FIRST  HEAD  GFAR  CONFIGURATION  ?“|Q3 
135  V3  “  03  »  03 
140  PRINT 

145  INPUT  “WHAT  IS  THE  MEAN  VALUE  OF  THE  FORWARD  ENDURANCE  I  I ME  FOR  THE  S 
ECOND  HEAD  GEAR  CONFIGURATION  ?"|A2 
148  PRINT 

150  INPUT  “WHAT  IS  1  HE  STANDARD  DEVIATION  OF  THE  FORWARD  ENDURANCE  TIME  F 
OR  THE  SECOND  HEADGEAR  CONFIGURATION  ?'|Q2 
•55  V2  -  Q2  *  02 
160  PRINT 

165  INPUT  "WHAT  IS  THE  MEAN  VALUE  OF  THE  LATERAL  ENDURANCE  TIME  FOR  THE  S 
ECOND  HELMET  CONFIGURATION  ?"(A4 
168  PRINT 

170  INPUT  "WHAT  IS  THE  STANDARD  DEVIATION  Of  THE  LATERAL  ENDURANCE  TIME  F 
OR  THE  SECOND  HEADGEAR  CONFIGURATION  2-,Q4 
175  V4  =»  U4  •  Q4 

300  II  •  VI  /  6:Z2  «  V2  /  6 

305  23  “  V3  /  6: Z4  =  V4  /  * 

310  W1  *>  (A1  -  A2>  /  SQR  ( Z 1  ♦  Z2*> 

315  W2  =  ADS  (Wl) 

317  W3  -  (A3  -  A4)  /  60K  (33  ♦  Z4> 

319  W4  «  ADS  (W3) 

350  PR*  1 
400  PRINT 
450  PRINT 

500  PRINT  “PUR  THE  FOLLOWING  COND I T 1 ONGi “ 

502  PRINT 

305  PRINT  “HEADGEAR  CONf  1 GURAT I  ON  »li“ 

810  PKINT  “FORWARD  ENDURANCE  TIME)  "lAlj"  SECONDS." 


D-l 


312  PRINT  * (BTD.  D*V.«  "|Q1|*  SECS. > “ 

313  PRINT  "HEADGEAR  CONP I»JRAT ION  62i" 

320  PRINT  FORWARD  ENDURANCE  T I HE i  •|A2|"  SEDONDS . 

322  PRINT  "  <STD.  DCV.=  ■■|Q2;"  SECS.)" 

525  PRINT 

330  PRINT  "THE  CALCULATED  WELCH  NUMBER!  “ | W2 

330  XI  ■  W2  -  2.01 

333  IF  XI  <  O  THEN  GOTO  600 

3*0  X 2  -  M2  -  2.37 

363  IF  X2  <  O  THEN  GOTO  610 

570  X3  »  W2  -  4.03 

575  IF  X3  O  THEN  GCTQ  £,20 

5S0  GOTO  630 

600  PRINT  "THIS  IS  NOT  SIGNIFICANT  < p  > . 05 ) . ” 

605  GOTO  700 

610  PRINT  "THIS  IS  SIGNIFICANT  AT  THE  957.  CONFIDENCE  LEVEL  <p<.05>." 

615  GOTO  700 

620  PRINT  "THIS  IS  SIGNIFICANT  AT  THE  97.57.  CONFIDENCE  LEVEL  lp<.025>." 
625  GOTO  700 

630  PRINT  "THIS  15  SIGNIFICANT  AT  THE  99.57.  CONFIDENCE  LEVEL  <p< ,005) . " 
635  GOTO  700 
700  PRINT 
702  PRINT 

704  PRINT  "FOR  THE  FOLLOWING  CONDITIONS!" 

706  PRINT 

706  FR1NT  "HEDGEAR  CONFIGURATION  61i" 

710  PRINT  "LATERAL  ENDURANCE  TIME!  "»A3l"  SECONDS. ' 

712  PRINT  " (STD.  DEV.-  "|03l"  SECS.)" 

715  PRINT  "HEADGEAR  CONFIGURATION  #2i  " 

720  PRINT  "LATERAL  ENDURANCE  TIMEi  " |A4|"  SECONDS." 

722  PRINT  "(STD.  DEV.«  "j04|"  SECS.)" 

725  PRINT 

727  PRINT  "THE  CALCULATED  WELCH  NUMBER I  " j W4 

730  X4  =  W4  -  2.01 

733  IT  X4  <  O  THEN  GOTO  800 

760  X5  -  W4  —  2.57 

765  IF  X5  <  O  THEN  GOTO  810 

770  X6  =  W4  -  4.03 

??M  IF  X6  0  THEN  GOTO  820 

780  GOTO  930 

800  PRINT  "THIS  IS  NOT  SlGG’^ TCANT  (p>.05)" 

803  GOTO  890 

810  PRINT  "THIS  IS  SIGNIFICANT  AT  THE  957.  CONFIDENCE  LEVEL  (p  .05)." 

813  GOTO  890 

820  PRINT  "THIS  IS  SIGNIFICANT  AT  THE  97.57.  CONFIDENCE  LEVEL  <p  .025>." 
025  GOTO  890 

030  PRINT  "THIS  IS  SIGNIFICANT  AT  THE  99.57.  CONFIDENCE  LEVEL  <p<.005>." 

835  GOTO  890 

870  PRINT 

093  PRINT 

879  PR#  O 

900  INPUT  "DO  YOU  WISH  TO  EVALUATE  ANOTHER  PAIR  OF  HEADGEAR  CONFIGURATION 
S?  TYPE  1'  FOR  YES'  OR  TYPE  2  FOR  NO  . " | P 1 
910  IF  FI  -  1  THEN  GOTO  50 
977  END 
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